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ABSTRACT 


The kinetics of evaporation in binary alloys have been quanti- 
tatively treated. The formalism so developed works surprisingly 
well for several systems so far studied, notwithstanding the sev 
eral major as sumptions involved. In this memorandum, we have 
studied the kinetics of purification of beryllium through evapora- 
tion data actually acquired during vacuum induction melting. This 
study shows that normal evaporation equations are generally valid 
and useful for understanding the kinetics of beryllium purification. 
The normal evaporation analysis has been extended to cover cases 
of limited liquid diffusion. It has been shown that under steady- 
state evaporation, the solute concentration near the surface may 
be up to six orders of magnitude different from the bulk concentra- 
tion. Corrections for limited liquid diffusion are definitely 
needed for the highly evaporative solute elements, such as Zn, Mg, 
and Na, for which the computed evaporation times are improved by 
five orders of magnitude. The commonly observed logarithmic rela- 
tion between evaporation time and final concentration further sup- 
ports the validity of our normal evaporation equations. 
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INTRODUCTION 

The kinetics of evaporation in binary alloys have been quanti- 
tatively treated by C. H. Li (Refs. 1,2), The formalism so devel- 
oped works surprisingly well for the several systems so far studied, 
notwithstanding the several major assumptions involved. Specifi- 
cally, in this normal evaporation approach, it is assumed that 

a) the evaporating alloy is always homogeneous in composition, 

b) the alloy follows Raoult's law, and c) the vapor is instantly 
removed. In this work, we have studied the kinetics of puri- 
fication of beryllium through evaporation data actually acquired 
during vacuum induction melting by R. F. Bunshah and R. S. Juntz 
(Ref. 3). 


1 



EXPERIMENTAL DATA 

The details of the vacuum induction melting procedure are as 

follows. Beryllium (SR grade flakes) in the form of a cylindrical 

compact was crucible -free induction melted under an ambient pres - 
"”6 

sure of 10 torr. The cylindrical compatct was l - 3/4 inches in 
diameter and 4 inches high, and weighed about 100-130 g. Since 
the actual temperature of the melt was not known, for the purpose 
of our computations, the temperature was assumed, for the sake of 
comparison, to be 1500 °C, a temperature commonly used in Be 
melting (Ref. 3) . 
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EQUATION OF NORMAL EVAPORATION 


The solutions of the kinetic equations for two important cases 
are given in papers by Li and Mukherjee et al. (Refs. 2,4). In 
Case I, the solute evaporating rate, U, is much greater than the 
solvent evaporating rate, V, i.e., U » V. The equation re~ 
lating the evaporating time t^ to final concentration, m> is, 
for this case. 


In 


(1 - m Jm 


m (1 

o 


m) 


m ” m o 

(1 " m Q ) (1 " m ) 


( 1 ) 


and in Case II when U « V, the evaporation time is given as 
follows : 


where 


t 2 G 2 


In 


(1 - m )m 
cr 

m (1 * m) 
o 7 


m - m 


nun 


( 2 ) 


G x = " N q (1 - m Q )/AU 

G 0 = + N m /AV, and 
2 o o 

m,m = mole fractions of solute at t = t and t = 0, 
respectively. 


G’s are the time constants of evaporation. For Case I when 
U » V, solute depletion takes place (or m < m Q ) and, therefore, 
the logarithmic term in the parentheses is negative (the second 
term in the parentheses is, for small values of m and m^, neg- 
ligible by comparison with the logarithmic term) . Since the evap- 
orating time t^ must be positive, G-^ must be negative. G 2 > on 
the other hand, is always positive. 
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to reach the final concentrations 


In Table 1 the times t . , 

x 

m from initial concentration m Q (through evaporation from area A 
and with an initial total of moles of both solute and solvent) 

for each element in the melt have been compiled using Eqs. (1) and 
(2) above. In the table, the ideal evaporation rates at temperature 
T°K of each element have been computed using the formula given in 
Ref. 5, viz . , 

for the solute: 

A ~B /T 

U = K 10 U U X /ITT 

v u 

and (3) 

for the solvent: 

A -B /T 

v = k io v v / yrnf , 

where K = 5.833 x 10 ^ for metals; A , B , A , and B are 

’ u u v v 

elemental evaporation constants given in Ref. 5; and the M’s are 
molecular weights. 


i 
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TYPES OF IMPURITIES 

Three types of impurities can be distinguished: 

Type I. - The impurities that have much higher evaporating 
rates U than the solvent Be. For these impurities, such as Zn, 

Na, and Mg, evaporation Eq. (1) should be used. 

Type II. - The impurities that have much lower evaporating 
rates U than the solvent Be. For these impurities, such as Fe, 

Cr, Ni, Cu, Ti, and Si, Eq. (2) should be used. 

Type III. - The impurity that evaporates at about the same 
rate as the solvent Be, i.e., U ~ V. Only a single element, 

Al, is in this group. For this impurity, a different evaporation 
equation in infinite series form should be used (Ref. 2) . 
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PRELIMINARY DATA ANALYSIS 


Since the exact temperature of the melt was not known, the 
exact evaporation rates could not be computed. However, since all 
the solute and solvent elements were evaporating from the same 
liquid'gas interface of a fixed area A for the same length of time, 
we can compute, from the measured m o ’s and m' s, the values of 
P, which is the product of the evaporating time and evaporation 
rate for a given solute with the help of Eqs . (1) and (2). The 
value of P should be constant for all the solutes. 


i _ ideal 
1 C 1 


-N (1 - m ) 
o cr 

AU. 

l 


In 


(1 _ m )m m - m 

Qj , o 

m o (l - m) (1 - m) (1 - m Q ) ' 


(4a) 


or 


. ,• “N (1 ■ m ) 

thj = t^ d 1 U. = — 2—7 2- 

111 i A 


In 


(1 - m ) 
J o 

m o (l -m) 


m 


m 


(1 - m) (1 - m Q ) 


= P , (4) 


where t^ is the estimated evaporation time computed separately 
for each element at the assumed temperature of 1500°C. Similarly, 
from Eq. (2) and the ideal evaporation rate of the solvent at 

the same temperature, we have: 


ideal 

t 2 


V ideal 


= P . 


In Table 1 are listed the value of P’s for each element. P 
is defined, as shown, such that 


P = t_^ x IL when < tL 


P = t. x V, when . 
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It can be seen from Table 1 that P is almost constant for the low 
vapor pressure (Type II) impurities except for copper. 

Because of our three assumptions of normal evaporation, and 
because of the assumed evaporating temperature of 1500 °C, the 
computed values of t ^ vary from element to element, even though 
they should be identically the same. The variations in t>, how" 
ever, are not very large so that the equations presented here are 
still useful. 

The high vapor pressure (Type I) impurities (Mg, Na, Zn) 
also have uniform P values, which are roughly one order of magni" 
tude higher than those for the low vapor pressure (Type II) ele- 
ments. The negative value of P for aluminum shows that aluminum 
cannot be treated by Eqs . (1) and (2) because for this element 
U ~ V and a different evaporation equation in infinite series form 
should be used. 
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STEADY- STATE EVAPORATION 


To explain the large variation in P values among different 
types of impurity elements, one should consider limited liquid dif- 
fusion. The concept of effective evaporation rates, U e and 
upon reaching a steady-state condition is now introduced. 

P - t V = constant (la) 

e e ' 

or 

t^V^ = P = constant = t g V e , (2a) 

where the subscripts i refer to ideal or normal evaporation val- 
ues, while the subscripts e refer to effective or nonnormal val- 
ues . 

In the following paragraphs, the concept of U g will be de- 
veloped. 

In normal evaporation, we assume that the evaporating alloy 

is always homogeneous, i.e., the liquid diffusion constant = °o. 

But because of the large (orders of magnitude) differences in 

evaporation rates of elements, the surface concentration of the 

solute m in the melt is not the same as the bulk concentration 
s 

m Q due to the limited liquid diffusion. Figure 1 depicts the 
situation (at t = 0) (i.e., initial), and after reaching steady 

states, for the case U » V. 

For U » V, the evaporating surface region soon becomes de- 
pleted of solute, and the surface concentration of the solute, m , 
will be different from the bulk concentration, m , specifically, 

m « m . After the initial transients, a steady-state for the 
s o J 

alloy system and evaporating conditions will be reached. The steady 
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Fig. la Initial State, 
t = 0 


Fig. lb Steady State, 
t * t 

s 


state concentration profile is shown in Fig. lb. The equation of 
the profile, i.e. , m versus x, can be derived according to a 
procedure used by Tiller et al. (Ref. 6) in their study of freez- 
ing with limited liquid diffusion, i.e.. 



( 5 ) 


where the receding rate R of the evaporating surface can be com- 
puted from U. , V., m , M , M , and material densities. 

-L o U. v 

In steady-state. 


and 


m U. 
s 1 


m U 
o e 




(1 - m g )V. = (I - m Q )U e , 
or 


(1 


m U. 

S 3- 


- m )V. 
s' i 



(6a) 


(6b) 


(6e) 
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Solving for m g . 


Since 


m 


m V. 

o x 


s (1 " m ) U. + m V. 
v o x ox 


(7) 


m U. 
s x 


m U , 
o e 


( 8 ) 


where U is the effective evaporation rate of the solute due to 
the change in surface concentration. Substituting the m from 

3 

Eq. (7) in Eq. (8) gives 


m U.V. 

U = — — n = — — 1 

e m u i (1 “ m )U. + m V. 
o O l ox 

Similarly, for the case V » U, 

(1 - m )V. = (1 - m )V 
s' x o' e 

or 


(9) 


( 10 ) 


1 - m V.U. 

XT — — V — 1 1 

e 1 - m i (1 - m )U. + m V 


= U 


o x 


( 11 ) 


Thus, it can be seen mathematically that V = U . This means 

that in both cases, U » V, or V « U, the effective value of 

U or V is the same. This equality means that in the steady-state, 

the surface concentration m does not change with time, as ex- 

s 

pected. 

Now consider Eq. (7) in some detail. Several cases can be 
distinguished; these are: 
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U » V and m ~ m ~ 0 
o — — 


m V. 

_ _ o l 

s (1 - m )U. 
v o' 1 


m V. 
_s = _i 

m U. 

o 1 


m o v i » U 1 or v7 « ra o 


m = 1 . 
s 


mV. « U. or — » m 
o i l o 


mV. mV. 

m _ Qi o i 

s (1 - m )U. U. 

v o' 1 1 


m ~ 0 . 
o — 
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c) 


then 


when 


mV, ~ U . or 

o t — r 



*v in • 

— o 



m ~ 0 . 
o “ 


(iid) 


The effective evaporation rate of the solutes, U , can be 
computed by Eq. (9) . The effective evaporation rates have been 
variably raised or lowered to equal the ideal evaporation rate of 
the solvent (Be) . In the case of high vapor pressure impurities 
like Mg, Zn, Na, this means a lowering of evaporating rates by 
about five orders of magnitude . 

The surface concentration, m , and the ratio of surface to 

s 

bulk concentration have been computed in Table 1 by Eq. (7) . These 
ratios must also equal . Now we have 


P - t U 
e e 


Vi 


The effective time t , to reach the final concentrations m 

e 

have been computed for all the impurities in Table 1; t g should 
be a constant. It can be seen that for the low vapor pressure im- 
purities, t is almost constant. Again, the t for Cu is one 
order of magnitude less than others in the low vapor pressure im- 
purity group. Among the high vapor pressure impurities, the t 
' " ' e 

is almost constant. It can be seen that t for the high vapor 
pressure impurities is one order of magnitude larger than that of 
low vapor pressure impurities. 
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NONIDEAL SOLUTIONS 


The normal evaporation equations can be further modified by 
introducing the concept of solvent “Solute interactions . If one 
considers deviation from the Raoul t' s law, it can be shown that the 
evaporation rates of solvent and solute will be changed by a new 
factor, the activity coefficient. For very dilute solutions, 
which we are considering, evaporation rates for solute and solvent 
can be given by U / = y^U^ and V* = V g , respectively, where y^ 
is the activity coefficient of the solute (Henry constant, in this 
case) . 

In general, such consideration should not affect the estimated 
time, t g , for the type I impurities (Fe, Ni, Cr, Ti, etc.), ex- 
cept in the case of copper. From the available phase diagram, it 
is evident that Be and Cu atoms have significant interactions 
with each other and form many compounds (Refs. 7,8). Thus the 
excess free energy for the system should be negative, resulting in 
a lower value of y p compared to the ideal solution case 
(^y^ u = 1.0) , In such a situation, is not very much higher 

than U' , and this leads to a greater effective time, t 

wU 6 y LU 

to reach the final concentration. The exact value of y_ is 

Cu 

difficult to evaluate because of the lack of data in this composi- 
tion range . 

In the cases of Zn and Mg, the solubility in beryllium is 
negligible (Ref. 8) and so y^ and y^ should be greater than 
one. This would lower the estimated time, t £ , meaning a better 
agreement with the other elements . 

It is important to note that exact calculation for the devia- 
tion from the ideal solution case is difficult because of the lack 
of data, and such estimations are purely qualitative. 
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CONCLUSIONS 


1. Normal evaporation Eqs. (1) and (2) are generally valid 
and useful for the study of beryllium purification. This agrees 
with a previous study on the results of Ni"Fe and Ni"Cr sys” 
terns (Ref. 4). 

2. The computed values of P 85 t.U. or t.V. are essen- 

li l i 

tially constant for the same type of solute elements, although dif- 
fering by one order of magnitude between different types of ele- 
ments . 

3. The normal evaporation analysis has been extended to cover 

cases of limited liquid diffusion- It is shown that under steady- 

state conditions, the solute concentration near the surface, i.e., 

m , may be up to five orders of magnitude different from the bulk 
s 

concentration, m . 

o 

4. Corrections for limited liquid diffusion are definitely 

needed for the highly evaporative solute elements Mg, Zn, and 

Na, for which the computed evaporation times, t , are improved 

© 

by three or four orders of magnitude. 

5. After these corrections, the computed times are comparable, 
even between different types of solute elements. 

6. The agreement among the computed times can be further im- 
proved by considering the deviation from the ideal solution condi- 
tion. Because of lack of data, it is difficult to give any quanti- 
tative calculation. 

7. Estimations of solute evaporation behavior are now well 
within one order of magnitude. 



8. The commonly observed logarithmic relation between evapo - 
ration time t and final concentration m further supports the 
validity of our normal evaporation equations. Such relations are, 
however, true only for widely different solute and solvent evapo" 
rating rates and for dilute solutions, as can be seen from Eqs. (1) 
and (2) . 
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